The Pfaffian solution of a dimer-monomer problem: 
Single monomer on the boundary 
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We consider the dimer-monomer problem for the rectangular lattice. By mapping the problem 
into one of close-packed dimers on an extended lattice, we rederive the Tzeng-Wu solution for a 
single monomer on the boundary by evaluating a Pfaffian. We also clarify the mathematical content 
of the Tzeng-Wu solution by identifying it as the product of the nonzero eigenvalues of the Kasteleyn 
matrix. 
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I. INTRODUCTION 

An outstanding unsolved problem in lattice statistics 
is the dimer-monomer problem. While it is known [| 
that the dimer-monomer system does not exhibit a phase 
transition, there have been only limited closed-form re- 
sults. The case of close-packed dimers on planar lat- 
tices has been solved by Kasteleyn Q and by Temperley 
and Fisher 0, 0], and the solution has been extended 
to nonorientable surfaces 0, 0] . But the general dimer- 
monomer problem has proven to be intractable 0. 

In 1974 Temperley @ pointed out a bijection between 
configurations of a single monomer on the boundary of a 
planar lattice and spanning trees on a related lattice. The 
bijection was used in Q to explain why enumerations of 
close-packed dimers and spanning trees on square lattices 
yield the same Catalan constant. More recently Tzeng 
and Wu made further use of the Temperley bijection 
to obtain the closed-form generating function for a single 
monomer on the boundary. The derivation is however in- 
direct since it makes use of the Temperley bijection which 
obscures the underlining mathematics of the closed-form 
solution. 

Motivated by the Tzeng-Wu result, there has been re- 
newed interest in the general dimer-monomer problem. 
In a series of papers Kong [Tol fTTL IT^ | has studied numer- 
ical enumerations of such configurations on m x n rectan- 
gular lattices for varying to, n, and extracted finite-size 
correction terms for the single-monomer [lfj and general 
monomer-dimer [TillT^] problems. Of particular interest 
is the finding [l(J that in the case of a single monomer 
the enumeration exhibits a regular pattern similar to that 
found in the Kasteleyn solution of close-packed dimers. 
This suggests that the general single-monomer problem 
might be soluble. 

As a first step toward finding that solution it is nec- 
essary to have an alternate and direct derivation of the 
Tzeng-Wu solution without the recourse of the Temper- 
ley bijection. Here we present such a derivation. Our new 
approach points to the way of possible extension toward 
the general single-monomer problem. It also identifies 
that, besides an overall constant, the Tzeng-Wu solution 
is given by the square root of the product of the nonzero 



eigenvalues of the Kasteleyn matrix, and thus clarifies its 
underlining mathematics. 



II. THE SINGLE MONOMER PROBLEM 

Consider a rectangular lattice C consisting of an array 
of M rows and N columns, where both M and N are odd. 
The lattice consists of two sublattices A and B. Since 
the total number of sites MN is odd, the four corner 
sites belong to the same sublattice, say, A and there are 
one more A than B sites. The lattice can therefore be 
completely covered by dimers if one A site is left open. 
The open A site can be regarded as a monomer. 

Assign non-negative weights x and y respectively to 
horizontal and vertical dimers. When the monomer is on 
the boundary, Tzeng and Wu 9] obtained the following 
closed-form expression for the generating function, 



M — l N - 1 



n n 



4a; 2 cos 2 



M+ 1 



Ay z cos 



AT+1 



(1) 



This result is independent of the location of the monomer 
provided that it is on the boundary. 

We rederive the result Q using a formulation which 
is applicable to any dimer-monomer problem. We first 
expand C into an extended lattice CJ constructed by con- 
necting each site occupied by a monomer to a new added 
site, and then consider close-packed dimers on Since 
the newly added sites are all of degree 1, all edges orig- 
inating from the new sites must be covered by dimers. 
Consequently, the dimer-monomer problem on £ (with 
fixed monomer sites) is mapped to a close-packed dimcr 
problem on which can be treated by standard means. 

We use the Kasteleyn method [2) to treat the latter 
problem. Returning to the single-monomer problem let 
the boundary monomer be at site sq — 0-,n) as demon- 
strated in Fig. la. The site so is connected to a new 
site s' by an edge with weight 1 as shown in Fig. lb. To 
enumerate close-packed dimers on C using the Kasteleyn 
approach, we need to orient, and associate phase factors 
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FIG. 1: (a) A dimer-monomer configuration on a 5 x 5 lattice 
£ with a single monomer at so = (1, 3). (b) The extended lat- 
tice £ with edge orientation and a phase factor i to horizontal 
edges, (c) The reference dimer configuration Co. 



to, edges so that all terms in the resulting Pfaffian yield 
the same sign. 

A convenient choice of orientation and assignment of 
phase factors is the one suggested by T. T. Wu [l3| . 
While Wu considered the case MN — even, the con- 
sideration can be extended to the present case. Orient 
all horizontal (resp. vertical) edges in the same direction 
and the new edge from s' to sq, and introduce a phase 
factor i to all horizontal edges as shown in Fig. lb. Then 
all terms in the Pfafhan assume the same sign. To prove 
this assertion it suffices to show that a typical term in the 
Pfaffian associated with a dimer configuration C has the 
same sign as the term associated with a reference config- 
uration Co • For Co we choose the configuration shown in 
Fig. lc, in which horizontal dimers are placed in the first 
row with vertical dimers covering the rest of the lattice. 
Then C and Co assume the same sign. 

The simplest way to verify the last statement is to 
start from a configuration in which every heavy edge in 
Co shown in Fig. lc is occupied by two dimers, and view 
each of the doubly occupied dimers as a polygon of two 
edges. Then the 'transposition polygon' (Cf. |2() formed 
by superimposing any C and Co can always be generated 
by deforming some of the doubly occupied edges into big- 
ger polygons, a process which does not alter the overall 
sign. It follows that C and Co have the same sign for any 
C. This completes the proof. 

Here we have implicitly made use of the fact that the 
monomer is on the boundary. If the monomer resides in 
the interior of C, then there exist transposition polygons 
encircling the monomer site which may not necessarily 
carry the correct sign. The Pfaffian, while can still be 
evaluated, does not yield the dimer-monomer generating 
function. We shall consider this general single-monomer 
problem subsequently [Tlj . 

With the edge orientation and phase factors in place, 
the dimer generating function G is obtained by evaluating 
a Pfaffian 



where A' is the antisymmetric Kasteleyn matrix of di- 
mension (MN + 1) x (MN + 1). Explicitly, it reads 
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Here, A is the Kasteleyn matrix of dimension MN for C 
given by 



A = ixT M <8> In + ylm ® T N , 



(4) 



with In the N x N identity matrix and T/v the N x N 
matrix 
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Note that elements of A are labeled by {(to, n); (ml , n')}, 
where (to, n) is a site index, and the clement 1 in the first 
row of A' is at position (1, n) of A, n = odd. 

Expanding along the first row and column, we ob- 
tain 

DetA' = C(A;{(l,n);(l,n)}) (6) 

where C(A; {(l,n); (1, n)}) is the cofactor of the 
{(1, n); (1, n)}-th element of Det A. 

The cofactor C(a,(3) of the (a, j3)-th element of any 
non-singular A can be computed using the identity 

C(A;a,/3) = A~ l {p,a) x DetA, (7) 

where A -1 ((3, a) is the (/?, a)-th element of A. However, 
the formula is not directly useful in the present case since 
the matrix A is singular. We shall return to its evaluation 
in Sec. IV. 

III. EIGENVALUES OF THE DETERMINANT A 

In this section we enumerate the eigenvalues of A. 

The matrix Jjy is diagonalized by the similarity trans- 
formation 



Uk 1 Tjv Un = A 
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where Un and U N are N x N matrices with elements 
U N (ni,n 2 ) = 
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and Ajv is an N x N diagonal matrix whose diagonal 
elements are the eigenvalues of Tjy, 



A™ = 2i cos ■ 



N+l' 



m = 1,2, • • • ,N. 



(9) 



Similarly the MN x MN matrix A is diagonalized by 
the similarity transformation generated by Umn = Um® 
Un- Namely, 



Umn A Umn — ^-mn 
where Amn is a diagonal matrix with eigenvalues 
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m = 1,2, ...,M, n = 1,2, ...,N. 
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on the diagonal, and elements of Umn and U M l N are 

UMN{mi 1 n 1 ;m 2 ,n 2 ) = U M (mi, ma) U N {n\, n 2 ) 
U M 1 N (mi,ni- 1 m2,n 2 ) = U^(m\, m 2 ) f7^ x (ni, n 2 ). 



Then we have 
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As in (J2J, close-packed dimers on £ are enumerated 
by evaluating VDetA For MN even, this procedure 
gives precisely the Kasteleyn solution 0. For MJV odd, 
the case we are considering, the eigenvalue l mn = for 
to = (M + l)/2,n = (AT + l)/2, and hence Det A = 0, 
indicating correctly there is no dimer covering of C How- 
ever, it is useful for later purposes to consider the product 
of the nonzero eigenvalues of A, 
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where the prime over the product denotes the restriction 
Using the identity 
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where the factor Q is the product of factors with either 
m = M±L or n = ^±±. Namely, 
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where we have made use of the identity 

jV-l 

un \ N+l 



n 4 



: COS 



A + 1 



, AT = odd. 



The expression 11141 for P will be used in the next section. 



IV. EVALUATION OF THE COFACTOR 

We now return to the evaluation of the cofactor 
C(A; {(1, n); (1, n)}). We shall however evaluate the co- 
factor C(A\ {(to, n); (to', n')}) for general to, to', n, n', al- 
though only the result of to = m' = 1, n — n' is needed 
here. 

To circumvent the problem of using J7J) caused by the 
vanishing of Det A = 0, we replace A by the matrix 

A(e) = A + ehiN , e^O 

whose inverse exists, and take the e — * limit to rewrite 
as 



C(A;{(m,n);(m',n')}) 



lim 



[A^ 1 (e)] (to', n'; to, n) x Det A(e) 
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Quantities on the r.h.s. of 11()|) are now well-defined and 
the cofactor can be evaluated accordingly. The consider- 
ation of the inverse of a singular matrix along this line is 
known in mathematics literature as finding the pseudo- 
inverse 0, ^| . The method of taking the small e limit 
used here has previously been used successfully in the 
analyses of resistance ^3 an< l impedance 01 networks. 

The eigenvalues of A(e) are l mn {e) = / m „ + e and hence 
we have 

M JV 

Det A{e) = JJ J] [l mn + e ]=eP + (9(e 2 ), (17) 

m—l n—1 



where P is the product of nonzero eigenvalues given by 

d- 

We next evaluate A~ 1 (e)(m.n; to, n) and retain only 
terms of the order of 1/e. Taking the inverse of l|lUfl with 
A(e) in place of A, we obtain 

A-\e) = Umn A M 1 JV (e) U^ N . 

Writing out its matrix elements explicitly, we have 
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For e small the leading term comes from (m",n") = ( M 2 ^ 1 , ^-f-) for which l m ",n" — 0. Using U MN (m,n;m' ,n') 
U^ I 1 (m,m / )U^ 1 (n,n') and (JSJ, this leads to the expression 
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Thus, after making use of 1|16(1 and Ijl7jl we obtain 
C(A;{(m,n);(m',n')}) 
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Finally, specializing to m — m! — 1, n = n' and com- 
bining ©, ©, and {EH), w e obtain 



G(x,y) = v/C(A;{(l,n);(l,n)}) 
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for n odd (A site) 



y (m + i)(jv + i) ' 

= 0, for n even (B site) 

(20) 

This gives the result after introducing (|14|l for P. It 
also says that there is no dimer covering if the monomer 
is on a B site. 

The expression (|20ll clarifies the underlining mathe- 
matical content of the Tzeng-Wu solution by iden- 
tifying it as the product of the nonzero eigenvalues of 
the Kasteleyn matrix. This is compared to the Kaste- 
leyn result 0] that for MN — even the dimer generating 
function is given by the product of all eigenvalues. 



V. DISCUSSIONS 

We have used a direct approach to derive the closed- 
form expression of the dimer-monomer generating func- 
tion for the rectangular lattice with a single monomer on 
the boundary. Our approach is to first convert the prob- 
lem into one of close-packed dimers without monomers, 
and consider the latter problem using established means. 
This approach suggests a possible route toward analyzing 
the general dimer-monomer problem. 

We have also established that the Tzeng-Wu solution 
is given by the product of the nonzero eigenvalues of 
the Kasteleyn matrix of the lattice. This is reminiscent to 
the well-known result in algebraic graph theory 19] that 
spanning trees on a graph are enumerated by evaluating 
the product of the nonzero eigenvalues of its tree matrix. 
The method of evaluating cofactors of a singular matrix 
as indicated by (|16|l . when applied to the tree matrix of 
spanning trees details of which can be easily worked out, 
offers a simple and direct proof of the fact that all cofac- 
tors of a tree matrix are equal and equal to the product of 
its nonzero eigenvalues. The intriguing similarity of the 
results suggests there might be something deeper lurking 
behind our analysis. 

I am grateful to W. T. Lu for help in the preparation 
of the manuscript. 
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